A PRIME ORBIT THEOREM FOR SELF-SIMILAR FLOWS AND 
DIOPHANTINE APPROXIMATION 
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Abstract. Assuming some regularity of the dynamical zeta function, we establish 
an explicit formula with an error term for the prime orbit counting function of a 
suspended flow. We define the subclass of self-similar flows, for which we give an 
extensive analysis of the error term in the corresponding prime orbit theorem. 



1. Introduction 

In IPr'Plll , Parry and Pollicott obtain a Prime Orbit Theorem for certain dynamical 
systems — the so-called 'suspension flows'. (See also [[PP2| , Chapter 6].) The first 
results of this kind were obtained in special cases by Huber |[Hu|| , Sinai |S|], and 



Margulis [^], among others. See ||PP1| , 2] and the relevant references therein, as 



well as the historical note in ||BKS|, p. 154]. Parry and Pollicott derive the first term 



in the asymptotic expansion of the counting function of prime orbits, by applying the 
Wiener-Ikehara Tauberian Theorem to the logarithmic derivative of the dynamical 
zeta function. An alternate approach was taken by Lalley in |[Lall| , 2], who considers. 



in particular, the (approximately) self-similar case. Using a nonlinear extension of 
the Renewal Theorem, he shows that in the nonlattice case, the leading asymptotics 
are nonoscillatory. In the lattice case, the leading asymptotics are periodic, and it 
becomes a natural question whether they are constant or nontrivially periodic. 
In a recent book |[LvF2|| , we have developed a theory of complex dimensions of 



fractal strings (one-dimensional drums with fractal boundary, see [[LP| , P^M|| ). These 



(geometric) complex dimensions — defined as the poles of the associated geometric 
zeta function — enable us to describe the oscillations intrinsic to the geometry or 
the spectrum of fractal drums, via suitable 'explicit formulas', obtained in |[LvI''2| , 
Chapter 4]. 

In this paper, we apply these explicit formulas to obtain an asymptotic expansion 
for the prime orbit counting function of suspension flows. The resulting formula in- 
volves a sum of oscillatory terms associated with the dynamical complex dimensions 
of the flow. We then focus on the special case of self-similar flows and deduce from 
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our explicit formulas a Prime Orbit Theorem with error term. In the lattice case 
(to be defined below), the counting function of the prime orbits, ip^ji^x), has oscilla- 
tory leading asymptotics and our explicit formula enables us to give a very precise 
expression for this function in terms of multiplicatively periodic functions. In the 
nonlattice case (which is the generic case), the leading term is nonoscillatory and we 
provide a detailed analysis of the error term. The precise order of the error term de- 
pends on the 'dimension free' region of the dynamical zeta function, as in the classical 
Prime Number Theorem. This region in turn depends on properties of Diophantine 
approximation of the weights of the fiow. 

For suspension fiows, the dynamical complex dimensions are defined as the poles 
of the logarithmic derivative of the dynamical zeta function. On the other hand, the 
geometric complex dimensions of a fractal string are defined in |[LvFl| , 2] as the poles 
of the geometric zeta function, which coincides with the dynamical zeta function when 
the string and the fiow are self-similar. Thus the geometric complex dimensions of 
a self-similar fiow only depend on the poles of the corresponding zeta function, and 
they are counted with a multiplicity, whereas the dynamical complex dimensions of 
a fiow depend on the zeros and the poles of the dynamical zeta function, and they 
usually have no multiplicity. Due to the fact that the dynamical zeta function of 
a self-similar fiow has no zeros, the two sets of complex dimensions coincide in this 
case. 



2. The Zeta Function of a Dynamical System 

Let > be an integer and let = {1, . . . , A^}^ be the space of sequences over 
the alphabet {1, . . . , N}. Let tr : — (0, oo] be a function, called the weight. On Q, 
we have the left shift a, given on a sequence (a„) by {aa)n = fln+i- We define the 
suspended flow J^n, on the space [0, oo) x as the following dynamical system (time 
evolution, see |Pr'P2| , Chapter 6]): 



(2.1) ^„(t,a) 



(t,a) ifO<t<tt)(a) 
jF„(t — tt) (a), era) if t > tt)(a). 



(Note that J-'nit, a) may not be defined. However, it is always defined on periodic 
sequences.) This formalism is seemingly less general than the one introduced in |PP2| , 
Chapter 1]. However, defining tt)(a) = oo when the sequence a contains a prohibited 
word of length 2, and e~*°° = 0, allows us to deal with the general case. 
Given a finite sequence j: = ai, a2, . . . , a; of length / = /(j;), we let 

a = ai, a2, . . . ,ai, ai, 02, . . . , a;, . . . 

be the corresponding periodic sequence, and we define crj = 02, . . . , a/, ai. The total 
weight of the orbit of a on j: is 

(2.2) tt)t(j:) = ro (a) + ra((Ta) H hra(o-'"^a). 

We now define (see [^,0 and [|PFg, Chapter 5]): 
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Definition 2.1. The dynamical zeta function of JF^ is defined as 



(2.3) Us) = ^w{^^^ 



wfiere tlie sum extends over all finite sequences f of positive lengtli. 

For = 0, the alphabet is empty, and we interpret JF„ as the static flow on a 
point, and Cro('S) = 1- Further, for = 1, we have the dynamical system of a point 
moving around a circle of length nit{l) = rD(l, 1, . . . ); Cm{s) = (1 — e~*"'**^^))~^. 

We also introduce the logarithmic derivative 



(2.4) -^W-E^ 



For > 1, this series does not converge for s = 0. We assume that (|2.4| ) converges 
for some value of s > 0, and the abcissa of convergence of this series will be denoted 
by D, the dimension of J^tn-Q Clearly, D > 0. Then ( p.4| ) is absolutely convergent 
for Res > D. Moreover, as in [ [LvFl| , 2], we assume that there exists a function 



S: M — > M, called the screen, satisfying S{t) < D for every t G M, such that — Cro/Cn 
has a meromorphic extension to a neighborhood of the region 

(2.5) W = {s = a + it: a>S{t)}, 

called the window. In Section^, we will also assume that —('ro/Cm satisfies the growth 
conditions (Hi) and (H2), to be introduced in Section ^. We will then say that J-'^ 
satisfies (Hi) and (H2). 

Definition 2.2. The poles of — Crt3/Cn3('S) iii ^ are called the complex dimensions of 
the flow JF„. The set of complex dimensions of JF„ in W is denoted by VniW) or 
for short. 

The nonreal complex dimensions of a flow come in complex conjugate pairs uj, ZJ 
(provided that W is symmetric about the real axis). If has a meromorphic ex- 
tension to W as well, then the complex dimensions of JF„ are simple and they are 
located at the zeros and poles of C^, 

V^W) = {u;eW: (U^) = or 00}, 

and the residue at a complex dimension u (i.e., Tes{—Cn/ Cm] ^)) " OTd{(n; u), where 
OTd{(n,; u) = n is the order of Cm at lu: Cn}{s) = C{s — u;)" + 0((s — c<j)"+^). In 
general, the complex dimensions of JF„ in W are not simple, and the residues are 
not necessarily integers. By abuse of notation, we write ord(Cro;'^) = T^^^iCm/Cm]^) 
if Cm/ Cm has a meromorphic extension with a simple pole at uo, even if the residue is 
not an integer (and consequently. Cm is not analytic at u). 



^The dimension often coincides with the topological entropy of the flow; see [PP2, Chapter 5] 
and the references therein. 
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2.1. Periodic Orbits, Euler Product. A periodic sequence a in Q with period /, 
a = ai, . . . ,ai,ai, . . . ,ai, . . . , gives rise to the finite orbit {a, era, . . . , cr'~^a} of cr. It 
is clear that / is a multiple of the cardinality #{ci, era, . . . , (j'~^a} of this orbit. 

Definition 2.3. A finite sequence f is primitive if its length /(p) coincides with the 
length of the corresponding periodic orbit of a. 

We denote by a\Q the space of periodic orbits of a. Thus 

(2.6) (7\il = {{(x'^f : A; G N} : p is a finite sequence} . 

We reserve the letter p for elements of a\Q. So p will denote a periodic orbit of a, 
and we write ^^p for its length. The total weight of an orbit p is 

(2.7) m,(p) = 5^m(a). 

Theorem 2.4 (Euler sum). For lies > D, we have the following expression for the 
logarithmic derivative of Cxo '■ 

/-/ °° 

(2.8) -f{s)= J2 E"'*(P)^"^'""^'^' 

p£a\n k=l 

where p runs through all periodic orbits of JF„. 



Proof. We write the sum in ( p.4|) over the finite sequences p as a sum over the primitive 
sequences and repetitions of these. An orbit p of cr contains different primitive 
sequences of length #p, so we obtain 



j: ^ ' r:primitivc k=l 



-ksvotif) 



oo 



ksWtip) 



p&cr\Q k=l 

The theorem follows. □ 

Definition 2.5. The following function counts the periodic orbits and their multiples 
by their total weight: 

(2.9) Mx)= Yl 

kWt{p)<logx 

The function is the counterpart of ipi^) = J2p''<x^'^SP^ which counts prime 

powers with a weight \ogp; see Example p.6| . 

Corollary 2.6. We have the following relation between Cro/Cw O'l^'d ipn '■ 

/■I poo 

(2.10) -^{s) = / x-^#«(x), 
for Re s > -D. 
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The integral on the right-hand side of (|2.1CI|) is a Riemann-Stieltjes integral asso- 
ciated with the monotonia function ip^. 

Corollary 2.7 (Euler product). The function Cro(s) has the following expansion as 
a product: 



(2.11) Us) = n 



srot(p) ' 

peo-\n 

where p runs over all periodic orbits of The product converges for lies > D. 
Proof. In ( p.8| ), we sum over k to obtain 

^ pecT\n pecT\n 

The theorem then follows upon integrating and taking exponentials. □ 

In Section ^, we combine the above Euler product representation of —Cm/Cm with 
our explicit formulas of Section ^ to derive a Prime Orbit Theorem for primitive 
periodic orbits. 

Remark 2.8. We use ipw instead of the more direct counting function 

^wix) = 1- 

rot(p)<loga; 

However, setting 6'„(x) = Emt{p)<ioga' "'t(P)' so that ?/'„(x) = 9n{x) + 6'„(x^/^) + 

m{x)j , as X ^ oo, we find 



^«(^) = / ^ — 7de„{t) - 



\ogt ' logx ' Jq log^t t ' 
from which it is easy to derive the corresponding theorems for from those for tp^- 

2.2. Self-Similar Flows. A self-similar flow is best viewed as the following dynamics 
on the region of Figure |l|. A point x = xiN'-^ + X2N~'^ + ■ ■ ■ = .xiX2 ... on the unit 
interval moves vertically upward with unit speed until it reaches the graph, at which 
moment it jumps to {Nx} = Nx — [Nx] = .X2X3 . . . , the fractional part of Nx, and 
continues from there. In Figure 0, = 5, and the expansions of A, B and C in base 
5 are A = 17/124 = .032, B = 85/124 = .320, C = 53/124 = .203. 

Definition 2.9. A flow J-"^ is self-similar if N > 2 and the weight function tn 

depends only on the first letter of the sequence on which it is evaluated. We then 
put 

(2.13) Wj = rv{j,j,j,...), 
and 

(2.14) r(j:) = e-"'*^, r, = e""'^' = r{j,j,j, ...), 

for j = 1, . . . , A^. The numbers rj are called the scaling ratios of JF^. 
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OA C B 

Figure 1. A self-similar flow, = 5, with the orbit of 17/124. 

Note that < rj < 1. We will assume that the weights Wj = logr^^ are ordered 
in increasing order, < wi < W2 < ■ ■ ■ < wn, so that 1 > ri > r2 > ■ ■ ■ > vn > 0. 
When N = 2, the flow is called a Bernoulli flow. Such flows play an important role 
in ergodic theory (see ||BKS| , Chapters 2, 6 and 8]). 



Theorem 2.10. The dynamical zeta function associated with a self- similar flow has 
a meromorphic continuation to the whole complex plane, given by 

(2.15) CM ^ 

Its logarithmic derivative is given by 
(2.16) 



c 

-fis) 



E 



j=l ^f j 



The dimension D > of the flow is the unique real solution of the equation 1 = 
Z^j=i ' J- 

Proof. The sum over periodic sequences of fixed length / can be computed as follows: 

N N N 

E '■(f)' = EE---E<.-'-'-: 



ai=l a2=l 



ai=l 



[rl-\ h r 



NJ 
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Hence, for Res > D, the sum over all periodic sequences is equal to 

1=1 y:/(y)=Z 1=1 ^ j=l ^ 

The theorem follows upon exponentiation and analytic continuation. Since the func- 
tion 1 — 'Yl!j=i holomorphic, C,xo is meromorphic. □ 



Remark 2.11. Because of Theorem |2.1CI| , for a self-similar flow we can take the full 
complex plane for the window, = C; in that case, there is no screen. However, in 
applying our explicit formulas, we sometimes choose a screen to obtain information 
about the error of an approximation. 

Corollary 2.12. The set of complex dimensions T>^ = Pto(C) of the self- similar flow 
JF„ is the set of solutions of the equation 

N 

(2.17) E^T = 1' 

Moreover, the complex dimensions are simple {that is, the pole of —('^/(^nj at uj is 
simple). The residue at uj equals — ord(Cn);^)- 

2.2.1. Connection with Self-Similar Fractal Sets. Given an open interval I of length L, 
we construct a self-similar one-dimensional fractal set £. with scaling ratios ri, . . . ,r^. 
Such a set is called a fractal string (see [[LP| , [LM| , [LvFl| , 2]). The following construc- 



tion is reminiscent of the construction of the Cantor set. Let scaling factors 
ri, r2, . . . , Tat be given (A^ > 2), with 

1 > ri > r2 > . . . > Tat > 0. 

Assume that 

N 

(2.18) i?:=J]r, <1. 

i=i 

Subdivide / into intervals of length riL, . . . yr^L. The remaining piece of length 
(1 — i?)L is the first member of the string, denoted by /i, also called the first length 
in Remark p^.l4| below. Repeat this process with the remaining intervals, to obtain 



A^ new lengths I2, . . . , In+i in the next step, and A^'^ ^ new lengths in the k-th step. 
As a result, we obtain a self-similar string C consisting of intervals of length L{1 — 
R)rl^ . . . r^ {ki, . . . , /ca? G N), and a sequence h > h ^ h ^ ■ ■ ■ of positive numbers, 
called the lengths of the string. We let (c{s) = Yl^i ^jy ^^e geometric zeta function 
of £ (see F^^ , 2]). 

Theorem 2.13. Let C be a self-similar string, constructed as above with scaling 
ratios ri = e~^^, . . . , r^v = e~^^ . Then the geometric zeta function of this string has 
a meromorphic continuation to the whole complex plane, given by 

(2.19) (c{s) = (L(l - R))XM, for s G C. 

Here, L is the total length of C, and R is given by ( ^.18[ ). 
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This follows from Theorem |2.10| combined with ||Lvi''2| , Theorem 2.3, p. 25] 



Remark 2.14. For a self-similar string, the total length of C is also the length of the 
initial interval I in the above construction. We can always normalize a self-similar 
string in such a way that C,c = Cm (equivalently, that the first length of £ is 1), by 
choosing L{1 — R) = 1. This does not affect the complex dimensions of the string. 



Note that we need to assume that R = Y2f=i '"i < 1; which corresponds to a lower 
bound on the weights wj = — logrj. There is no analogue of this condition for general 
suspended flows. 

Remark 2.15. The Euler product does not seem to have a clear geometric interpre- 
tation in the language of fractal strings. There is, however, a natural self-similar flow 
on C: the flow 



(2.20) Mt,J,x) 



(0, j, xe*) if xe* < Ij, 

J-'c{t — log I j, j, 1) otherwise. 

The lengths Ij correspond to the periodic sequences j: of the flow JF„ via the formula 

«(y)-i 

(2.21) = n 



k=0 



Remark 2.16 (Geometric and dynamical complex dimensions). In |[LvF2|| , the geomet 



ric complex dimensions of a fractal string are defined as the poles of its geometric zeta 
function. Thus the complex dimensions are counted with a multiplicity, and the zeros 
of the geometric zeta function are unimportant. On the other hand, the dynamical 
complex dimensions are defined as the poles of the logarithmic derivative of the dy- 
namical zeta function. Thus the complex dimensions are simple, and both the zeros 
and the poles of the dynamical zeta function are counted. For self-similar flows, the 
dynamical zeta function and the geometric zeta function of the corresponding string 
coincide (up to normalization), and this zeta function has no zeros. Hence, as sets 
(without multiplicity), the geometric and dynamical complex dimensions coincide for 
self-similar flows and strings. 

Remark 2.17 (Higher-dimensional case). We have discussed above the case of fractal 
strings (i.e., the one-dimensional case) because it is the one studied in most detail 
F2|| . However, it is clear that our results can be applied to higher-dimensional 



m 



self-similar fractals |[F|, |Mn|| as well. This allows us to obtain information about the 



symbolic dynamics of self-similar fractals. On the other hand, as in the previous 
remark, it does not give information about the actual geometry of such fractals. 

2.3. The Lattice and Nonlattice Case. Let JF„ be a self-similar flow. Recall that 
to depends only on the first symbol and Wj = tt)(j, j, . . . ) for j = 1, . . . , iV. Consider 
the subgroup G of M generated by these weights, G = J2j=i '^Wj. 

Definition 2.18. The case when G is dense in R is called the nonlattice case. We 
then say that JF^ is a nonlattice flow. 
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The case when G is not dense (and hence discrete) in R is called the lattice case. 
We then say that JF„ is a lattice flow. In this situation there exists a unique positive 
real real number w, called the generator of the flow, and positive integers ki, . . . ,k]sf 
without common divisor, such that 1 < ki < ■ ■ ■ < kjsf and 

(2.22) Wj = kjW, 
forj = l,...,iV. 

The generator of JF„ generates the flow in the sense that the weight of every periodic 
sequence of JF„ is an integer multiple of w. 

We introduce a real number Dq as follows: Let m be the number of integers j in 
1, . . . , such that rj = r^, and let Dq G M be deflned by 

N-m 

(2.23) 1 + rf" = mr^\ 

3=1 

The dynamical complex dimensions of a self-similar flow are described in the following 
theorem. For brevity, we will usually refer to them as the complex dimensions of 

Theorem 2.19. Let JF„ he a self- similar flow of dimension D and with scaling ratios 
1 > ri > ■ ■ ■ > r^yf > 0. Then the value s = D is the only complex dimension of JF„ 
on the real line, all complex dimensions are simple, and the residue at a complex 
dimension {i.e., res(— C^/Cn,; cj)) is a positive integer. The set of complex dimensions 
in C {see Remark |2.11| ) of JF„ is contained in the bounded strip Dq <R.es < D : 

(2.24) V„ = V„{C) c{seC: Do<Res<D}. 

It is symmetric with respect to the real axis and infinite, with density bounded by 

(2.25) # (!)„ n {cu G C: I lmuj\ < T}) < —T + 0(1), 

71 

as T oo. 

In the lattice case, Cfo{s) ^-5 a rational function of e"'"*, where w is the generator of 
So, as a function of s, it is periodic with period 27Ti/w. The complex dimensions 
uj are obtained by finding the complex solutions z of the polynomial equation {of 
degree /cat) 

AT 

(2.26) z^^ = 1, wtth e"""^ = z. 

Hence there exist finitely many poles uji{= D),uj2, ■ ■ ■ , such that 

(2.27) Vxo = {uJu + 2mn/w. n E Z,u = 1, . . . ,q}. 

In other words, the poles lie periodically on finitely many vertical lines, and on each 
line they are separated by2'K /w. The residue of the complex dimensions corresponding 
to one value of z = e"""^ is the multiplicity of z as a solution of ( |2.26| ). 

In the nonlattice case, D is simple and is the unique pole of („, on the line 
lies = D. Further, there is an infinite sequence of complex dimensions of T^i coming 
arbitrarily close {from the left) to the line Hes = D. There exists a screen S to the 
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left of the line Res = D, such that — C^/C^ satisfies (Hi) and (H2) with k = 
{see Equations ( |3.2| ) and (|3.3| ) below), and the residue of —('^/(^ at the pole uj in 
W is equal to 1. Finally, the complex dimensions of JF„ can he approximated {via 
an explicit algorithm, as described in |[Lv.L''2| , §2.6]) hy the complex dimensions of a 



sequence of lattice strings, with smaller and smaller generator. Hence the complex 
dimensions of a nonlattice string have an almost periodic structure. 

Corollary 2.20. Every self-similar flow has infinitely many complex dimensions 
with positive real part. 



Proof of Theorem ^.19{ . For a proof of these facts, see |[Lvl^'2| , Theorem 2.13, pp. 37- 



40] . The density estimate (p.25|) follows from the fact that the right-hand side of (p.25| ) 
gives the asymptotic density of the number of poles of Cto) counted with multiplicity. 
The 0(1) estimate improves | |Lvl''2| , Theorem 2.22, p. 47]. It is proved in | |Lvl''3[ | . □ 

2.4. Examples of Complex Dimensions of Self-Similar Flows. 

Example 2.21 (The Cantor flow). This is the self-similar flow on the alphabet {0, 1}, 
with two equal weights Wi = W2 = log 3. It has 2" periodic sequences of weight n log 3, 
for n = 1,2,.... The dynamical zeta function of this flow is 

(2.28) Ccf(s) = 

After taking the logarithmic derivative, one finds that the dynamical complex dimen- 
sions are the solutions of the equation 

(2.29) 2-3~'^ = l (cjgC). 
We find 

f 27ii 

(2.30) Vcf = \D + n:neZ 



w 



with D = log3 2 and w = log 3. 



Example 2.22 (The Fibonacci flow). Next we consider a self-similar flow with two 
lines of complex dimensions. The Fibonacci flow is the flow Fib on the alphabet 
{0,1} with weights Wi = log 2, W2 = 2 log 2. Its periodic sequences have weight 
log 2, 2 log 2, . . . , n log 2, . . . , with multiplicity respectively 1,2,..., Fn+i, . . . , the Fi- 
bonacci numbers. Recall that these numbers are defined by the following recursive 
equation: 

(2.31) = Fn + with Fo = 0, F^ = 1. 

The dynamical zeta function of the Fibonacci flow is 

(2-32) CFib(^) = ^_2-]_^-s - 

The complex dimensions are found by solving the quadratic equation 

(2.33) {2-'^)'^ + 2-^ = 1 {ueC). 
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We find 2-^ = (-1 + ^5) /2 = 0-^ and 2"^ = -0, where = (1 + \^)/2 is the 
golden ratio. Hence 

D + n: n eZ\VJ {-D + (n + 1/2): n G Z 

w J I 

with D = log2 and w = log 2. 

Example 2.23 (The Golden fiow). We consider the nonlattice fiow GF with weights 
wi = log 2 and W2 = 01og2, where = (1 + \/5)/2 is the golden ratio. We call this 
fiow the golden flow. Its dynamical zeta function is 

(2-35) Cgf(.) = 

and its complex dimensions are the solutions of the transcendental equation 

(2.36) 2-" + 2-^" = l (ueC). 

A diagram of the complex dimensions of the golden fiow is given in Figure ^ To 
obtain it, we chose the approximation ^ 987/610 to approximate the fiow by the 
lattice fiow with weights wi = 610w, W2 = 9S7w, where w = (1/610) log2. We 
then used Maple to solve the corresponding polynomial equation. In particular, the 
dimension D of the golden fiow is approximately equal to D = .77921 .... See also 
Example |6.8| . 

3. Explicit Formulas 



We will formulate our explicit formulas in the more general framework of |P^vF2| , 
Chapter 4]. We refer to this book for the proofs and much additional information. 
Let ?7 be a positive measure on (0, oo), supported away from 0. Its Mellin transform 

is 

POO 

(3.1) Us) = / x-'drj, 

Jo 

the geometric zeta function of rj. We assume that Cr? is convergent for some s, and we 
write D for the abcissa of convergence. We assume that there exists a screen to the 
left of Re s = D such that C,ri has a meromorphic continuation to the corresponding 
window. To simplify the exposition, we also assume that the poles of Cr? are simple. 
This is the case, for example, if C,xv has a meromorphic continuation to a neighborhood 
of W . The general case, when the poles of C,ri niay have arbitrary multiplicities, is 
treated in |P^vF2| , Chapter 4]. 



The screen S is given as the graph of a bounded function S, with the horizontal 
and vertical axes interchanged: 

S={S{t)+it: t e R}. 

We will write inf S" = inftgK S'(t) and sup S* = sup^gjj S'(t). We assume in addition 
that is a Lipschitz continuous function; i.e., there exists a nonnegative real number, 
denoted by ||5'||Lip, such that 

\S{x)~S{y)\ < ||^||Lip|x-i/| forallx,i/GR. 
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Figure 2. The almost periodic behavior of the complex dimensions 
of the golden flow. 
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Further, recall from Section ^ that the window W is the part of the complex plane 
to the right of S; see formula (|2.5|). 

Assume that Ct? satisfies the following growth conditions: 

There exist real constants k > and C > and a sequence {Tn}nez of real 
numbers tending to ±oo as n ^ ±oo, with T_„ < < T„ for n > 1 and 
\imn-,^ooTn/\T^n\ = 1, such that 

(Hi): For all n G Z and all a > S(T„), 

(3.2) |C,(a + ^TO|<C■|TJ^ 
(H2): For all t G M, \t\ > 1, 

(3.3) \Cr,iS{t) + it)\<C-\t\\ 

Hypothesis (Hi) is a polynomial growth condition along horizontal lines (necessar- 
ily avoiding the poles of C^), while hypothesis (H2) is a polynomial growth condition 
along the vertical direction of the screen. 

In the following, we denote by Tes{g{s);uj) the residue of a meromorphic function 
g = g{s) at u. It vanishes unless is a pole of g. Also, 

(3.4) (.). = -j^, 

for k eIj. Thus, (s)o = 1 and, for A; > 1, (s)^ = s(s + 1) . . . (s + A; — 1). 
Let 

(3.5) N,{x) = ArW(x) = r]{0,x) + ^r/({x}), 

and more generally, let nI^''\x) be the {k — l)-st antiderivative of this function, for 
k = 1,2,.... Our explicit formula expresses this function as a sum over the poles of 

Theorem 3.1 (The pointwise explicit formula). Letr] be a generalized fractal string, 
satisfying hypotheses (Hi) and (H2). Let k be a positive integer such that k > k + 1, 
where k, > is the exponent occurring in the statement of (Hi) and (H2). Then, for 
all X > 0, the pointwise explicit formula is given by the following equality: 

N?{x)= Yl res (Us); u;)^^^ 

(^3^5) c^^{o,-i,...-(fc-i)} 



X) 



3=0 

Here, for x > 0, R{x) = -R^f' (x) is the error term, given by the absolutely convergent 
integral 

(3.7) «W = 4"W = i//'""«\t 
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Further, for all a; > 0, we have 

(3.8) = i?fl(x) < ||r||Lip)-^^-^maxK"P^,x^°f^} + C", 

' k — K, — 1 

where C is the positive constant occurring in (Hi) and (H2) and C is some suitable 
positive constant. The constants C{1 + ||r||Lip) and C depend only on t] and the 
screen, but not on k. 

In particular, we have the following pointwise error estimate: 

(3.9) R{x) = Rf{x) = O (x^'^P^+^-i) , 

as X ^ 00. Moreover, if S{t) < sup S for all t E M. {i.e., if the screen lies strictly 
to the left of the line Res = sup 5), then R{x) is of order less than ^'^^P'^+'^-i as 
X ^ 00: 

(3.10) R{x) = Rf{x) = o , 
as X 00. 

To formulate our second explicit formula, a distributional formula, we view r/ as a 
distribution, acting on a test function ip defined on (0, 00) by 

(3.11) {Nf\^)= / ^d^. 

Jo 

\k] 

We then define iV^ as the distribution obtained by integrating this one k times, so 
that 

It is easily verified that this definition coincides with formula ( |3.5| ) and the next line 
above when A; > 1. For k < 0, we extend this definition by differentiating \k\ times 



the distribution A'',^ . 



We shall denote by (p the Mellin transform of a (suitable) function if on (0, 00); it 
is defined by 



(3.13) ^(s) = / ip{x)x-'''dx. 

Jo 

Theorem 3.2 (The distributional explicit formula). Let rj be a generalized fractal 
string satisfying hypotheses (Hi) and (H2). Then, for every A; G Z, the distribu- 
tion nIj'^ is given by formula ( |3.6| ), interpreted as a distribution. That is, the action 
of nIj'^ on a test function tp is given by 

{Nf\^)= E res(a.);u;)^^^ 

(^3;L4^ c^^{0-l,...,-(fc-l)} 

fc-1 
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Here, the distribution R = Ri^ is the error term, given by 

(3.15) (R,^) = (i?W,v.> = ^JUs)^{s + k)^^. 

Definition 3.3. We will say that a distribution R on (0, oo) is of asymptotic order at 
most (respectively, less than x°) — and we will write R{x) = 0{x°') (respectively, 
R{x) = o(x")), as X — > oo — if applied to a test function (p, we have that 

(3.16) {R,(pa) = O (a°) (respectively, {R,ipa) = o(a")), as a ^ oo, 
where ipa{x) = a~^(p{x/a). 

Theorem 3.4. Fix G Z. Assume that the hypotheses of Theorem \3.3i are satisfied, 
and let the distribution R = be given by ( |3.15| ). Then R is of asymptotic order 
at most x™P'^+^~^ as X oo : 

(3.17) i?|fl(x) =0(x^"P^+'=~^), asx^oo, 

in the sense of Definition p. j| . 

Moreover, if S(t) < sup S for all t &W {i.e., if the screen lies strictly to the left of 
the line Re s = sup S), then R is of asymptotic order less than x™p^~^^~^ as x ^ oo: 

(3.18) R^\x) = a (x^"P^+'=-i) , as X ^ oo. 

We refer to [ [Lvl''2| , Chapter 4] for a proof of Theorems |3.1| , |3.2| and p.4| . 

Remark 3.5 (Oscillatory terms in the explicit formula). Our explicit formulas give ex- 
pansions of various functions associated with a measure as a sum over the poles 
of C»7- The term corresponding to the pole uj of multiplicity one is of the form Cx^ , 
where C is a constant depending on uj. If is real, the function x^ simply has a 
certain asymptotic behavior as x — > oo. If, on the other hand, = /9 + 27 has a 
nonzero imaginary part 7, then x^ = x^ ■ x*"^ is of order O(x^) as x 00, with a 
multiplicatively periodic behavior: The function x*"^ = exp(i7logx) takes the same 
value at the points e^'^^^'^'x (n E Z). Thus, the term corresponding to u will be called 
an oscillatory term. If there are poles with higher multiplicity, there will also be 
terms of the form Cx'^(logx)"*, m G N*, which have a similar oscillatory behavior. 

Example 3.6 (The classical Prime Number Theorem). Let ({s) = l+2"^+3^^+. . . 
(for Re s > 1) be the Riemann zeta function. This function has an Euler prod- 
uct C(s) = Ylpi^ ~ p~'')~^ (for Res > 1), where p runs over the prime numbers. 
Analogously to Corollary |2.6|, we obtain —('/Ci^) = x~^ dip^x) = (vi^)' "where 
V'(a^) = T.pk<x^ogp = N^^\x), and 



(3.19) 



m>l,p 
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is the prime string (see |[LvF2|| ). We apply Theorem |3.2| to rj = ^ to obtain the 
exphcit formula for ip: 

tP 1 ft' r/s 

(3.20) ^(x) =x-Y, MC/Cis); P)--7r- y{s)x'-, 

where p runs through the sequence of critical zeros of C,: C(p) = 0, < Rep < 1. 

By means of classical arguments |l|, Theorem 19], it is known that C has a zero free 
region of the form 

{a + it eC: a > 1-C/logt}, 

for some positive constant C. Also, — CVC is not too large on the boundary of this 
region. In our language, this means that we can choose a screen to the left of Re s = 1 
such that there are no zeros of C in and — C'/C satisfies (Hi) and (H2). Then (p.20|) 
becomes 

(3.21) il){x) = x + o{x), 

as X 00. This is equivalent to the classical Prime Number Theorem. 

Using the existence of the zero free region, one can derive the stronger estimate 



(3.22) i){x) =x + (^xe 



cyTogx 



as X —>■ 00, for some positive constant c (see [0; |, Theorem 23]). This is the classical 
Prime Number Theorem, with Error Term. 

4. The Prime Orbit Theorem for Flows 



Let J-'m be a suspended flow as in Section |^. In Corollary |2]^, we have written 
the logarithmic derivative of Cw('S) as the Mellin transform of the counting function 
ipw of the weighted periodic orbits of a, as defined in (pl9|). Put rj = dipro, so that 
Crj = — Cm/Cw- The poles of —(l^/Cm are the complex dimensions of JF„ and the 
residue at uj is —OTd{(n',^)- By Theorems |^ and we obtain the following 
explicit formula for the counting function of weighted periodic orbits of a. 

Theorem 4.1 (The Prime Orbit Theorem with Error Term). Let JF„ be a suspended 
flow that satisfies conditions (Hi) and (H2). Then we have the following equality be- 
tween distributions: 

(4.1) ^^(x) = ^+ J2 -ord(a;a;)- + resf-^^^;o')+i?(x), 

where ord ((^„; u) < denotes the order of at u, and 

(4.2) R{x) = - [ j^{s)x' — = (x™P^) , 



'S Cro ■S 
as X ^ 00. 

// is not a complex dimension of the flow, then the third term on the right-hand 
side of (|4.1|) simplifies to — Cro/Cro(0). In general, this term is of the form p + qlogx, 



for some constants p and q. 
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// D is the only complex dimension on the line lies = D, then the error term, 

is estimated by o{x^), as x —>■ oo. // this is the case, then we obtain a Prime Orbit 
Theorem for JF^ as follows: 

(4.4) ^^(a;) = ^ + o(x^), 

as a; oo. 

Proof. The first part of the theorem follows from the distributional explicit formula 



(Theorem p.2|) and from the first part of Theorem 3]^, while the second part follows 
from the second part of Theorem p.4|. □ 



The explicit formula holds for every flow satisfying our conditions on —('n/Cw- In 
particular, we can apply it to the 'axiom A flows' considered in |PP2| , Chapter 6], in 
view of [|PP2| , pp. 100-101]. We hope to do so more explicitly in a later work. For 
simplicity and for the sake of concision, however, we will focus in the rest of this 
paper on the important example of 'self-similar flows' (in the sense of Section |2l2 
above). 

5. The Prime Orbit Theorem for Self-Similar Flows 

For self-similar flows, does not have any zeros (see (p.3|) ). Hence every contri- 
bution to (|4.1| ) comes from a pole of (xo, and each coefficient — ord(CTO; t^) is positive. 
Furthermore, is never a complex dimension, so the third term on the right-hand 
side of (|4.1|) in the explicit formula is 



r 



We can obtain information about ip^v by choosing a suitable screen. 

5.1. Lattice Flows. In the lattice case, we obtain the Prime Orbit Theorem for 
lattice self-similar flows: 

1 ^ 

(5.2) = g,i\ogx)x'' - Y.'^. + O (a;^-") , 

i=i 

as X — s> oo. Here, D — a is the abcissa of the first vertical line of complex dimensions 
next to D, and the periodic function gi, of period w, is given by0 



(5.3) = „lo , = ^6r<»'"> 

D + Imn/w Oi — 1 



u -h zmn/w Oi — i ' 

n=— oo 



^We use the notation {u} for the fractional part, and [u] for the integer part of the real number 
u, so that {u} = u— [u] e [0, 1). 
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where bi = e^^. By choosing a screen located to the left of all the complex dimensions 
of ^ro, we can even obtain more precise information about ip^. In the notation of 
Theorem 12.191, we obtain 



i/^nix) = ^ - ord (Ctt,; uJu) ^ 



X 



(5.4) 



u=l 
Q 



UJu + 2TTin/w N 



1 ^ 



N 



= ^ - ord (Ctt,; UJu) ^„(loga;)x'^" - ^ ^ Wj 

where for each u 
(5.5) gu{y) = 



u=l j=l 

1, . . . , g, the function is periodic of period w, given by 

^2niny/w 



UJu + 211171/10 bu — I ^ ' 



where bu = e""^". Here, uji{= D),uj2, ■ ■ ■ ,ujq are given as in the lattice case of Theo- 
2.19| and ord(Ctn; 1^1) = ~1- 
For instance, the Cantor flow (with D = log^ 2 and Wi = W2 



rem 



w 



log 3, see 
Example ^.21j ) has 

(5.6) ^cf{x) = gi{logx)x'^ - 2 log 3, 

with gi{y) = 102^^^^^"^^ while the Fibonacci flowQ of Example p^.22| (with D = logg 
and Wi = w = log 2, W2 = 2w) has: 

(5.7) 7/'Fib(x) = ^7i(logx)a;^ + (72(logx)x"*/"x-^ - 3 log 2, 
where (?!(?/) = w(l)'^~^y^^^ and 



^2(y) = XI 



^2iTiny/w 



-D + 27ri(ri + l/2)/w 



In the second asymptotic term, the product e '^^{(loga:)/"'}^;'^*/'!" combines to give the 
sign (— l)[(iog^)/"']. 

5.2. Nonlattice Flows. In the nonlattice case, we use Theorem ^.19 according to 
which there exists 5 > and a screen 5* lying to the left of the vertical line Re s = D—5 
such that — C^o/Ctt) is bounded on 5* and all the complex dimensions uj to the right of 
S have residue res(— C^o/Cn,; 1^) = 1- Then R{x) = 0{x^~^), as x — > cxd. There are 
no complex dimensions with lie uj = D except for D itself. Hence, the assumptions 
of Theorem are satisfied. Therefore, in view of Theorem [4.1| , we deduce by a 
classical argument (see the proof of Theorems and |6.14| on page ^) the Prime 
Orbit Theorem for nonlattice suspended flows: 



D 



(5.8) 



ujev„\{D} 



n-s 



D 



D 



■^Also called the golden mean flow in the literature (see, e.g., [BK£, p. 59]), but not to be confused 
with the golden flow in our present paper, which is a nonlattice self-similar flow. 



A PRIME ORBIT THEOREM FOR SELF-SIMILAR FLOWS 



19 



as X oo. (See also Theorem p.l4| , and when N = 2, Theorem |6.7| below for a 



better estimate of the error.) We note that this estimate is always best possible, 



since by the nonlattice case of Theorem 2.19, there always exist complex dimensions 



of tt) arbitrarily close to the vertical line Hes = D. 

Remark 5.1. It would be interesting to apply Theorem to suspended flows that 
are more general than self-similar flows: for example, those considered by Lalley 
in ||Lall| , 2], such as the 'approximately self-similar flows' naturally associated with 
limit sets of suitable Kleinian groups. This would require a more detailed study of the 
dynamical zeta function of each of these flows. It is known that the lattice-nonlattice 
dichotomy applies in these more general cases; see ||Lall| , 2]. We hope to investigate 
this question in a later work. 

6. The Error Term in the Nonlattice Case 

A nonlattice flow has weights Wi < ■ ■ ■ < w^, where at least one ratio Wj/wk is 
irrational. Let 

N 

(6.1) /(,) = i_^e-'"^^ 

i=i 

Then D > is the unique real solution of the equation /(s) = 0. Moreover, the 
derivative 

N 

(6.2) f[s) = J2^,e--'' 

j=i 

does not vanish at D. 

When N = 2, the flow is called a Bernoulli flow (see Section |2.2| ). We write 
a = W2/W1, so a > 1 is irrational. In this case, we obtain very detailed information 
about the growth of —Cw/Cn on the line Res = D, and we can compute a pole free 
region for this function, by considering the continued fraction expansion of a. We 
briefly collect here the facts that we will use. See ||HaW| , P| , and ||Ba| , [vFl| , |vF2|| for a 
connection with the Riemann Hypothesis. 

6.1. Continued Fractions. Let a be an irrational real number with a continued 
fraction expansion a = [[oq, 01,02, . . .]] = + l/(ai + 1/(02 + •• • ))• We recall that 
the two sequences Oq, cti, . . . and ao, en, . . . are defined by ao = and, for n > 0, 
On = [«n], the integer part of a„, and On+i = l/(«n — ^n)- The convergents of a, 

(6.3) — = [[ao,ai,a2, . . . ,a„]], 

Qn 

are successively computed by p_2 = q-i = 0, p_i = g_.2 = 1, and Pn+i = ctn+iPn + 

Pn-i, Qn+i = cin+iqn + <in-i- We also define = cti ■ 0:2 and note the formula 

q'n+i = oin+iQn + Qn-i- Then 

(-1)" 

(6.4) g„a -p„ = — ; . 

1n+l 
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We have g„ > cp"' ^, where = (1 + v^)/2 is the golden ratio. 

Let n eN and choose / such that qi+i > n. We can successively compute (see |D|) 

n = diQi + ni, ni = di^iQi^i + rii^i, . . . , m = d^qo, 

where d^ is the quotient and n^, < is the remainder of the division of n^+i by q^- 
We set dij^i = di+2 = . . . = 0. Then 

oo 

(6.5) n = y^^d^q^. 

We call this the a-adic expansion of n. Note that Q < dy < a^j^i and that if d^ = a^^^i, 
then d^^i = 0. Also d^ < ai. It is not difficult to show that these properties uniquely 
determine the sequence dQ,di, . . . of a-adic digits of a. 

Lemma 6.1. Let n be given by ( |6.5|) . Let k > be such that dk ^ and dk-i = 
■ ■ ■ = do = 0. Put m = Y^^=kd'uPu- Then na — m lies strictly between {—lY I q'f,j^2 '^^^ 

Proof. We have na — m = Yl'^=k^'^(^1i' ~ P^)- Since aq^ — Pu = {—^T/qI+i: the 
terms in this sum are alternately positive and negative, and it follows that na — m 
lies between the sum of the odd terms and the sum of the even terms. To bound 
these terms, we use d^ < Ojy+i. Moreover, dk > I, hence dk+i < ak+2 — 1- It follows 
that na — m lies strictly between 

ak+i{aqk - Pk) + ak+^{,aqk+2 - Pk+2) + afe+5(a9fc+4 - Pk+i) + ■■■ 

and 

{aqk - Pk) + {ak+2 - l)(agfc+i - Pk+i) + 0^+4(05^+3 - Pfc+s) + 

Now Oy^i^aqy — Pu) = a^q^+i — qu^i) — {Pu+i —Pu-i)- So both sums are telescopic. 
The first sum immediately evaluates to —aqk-i +Pk-i = ll'k- The second sum 
equals {aqk - Pk) - {aqk+i - Pk+i) - {aqk - Pk) = i-'^)'' / Qk+2- ^ 

6.2. Two Generators: the Bernoulli Flow. Assume that N = 2, and let / be 

defined as in ( |6 . 1| ) with weights wi and W2 = awi, for some irrational number a > 1. 
We want to study the complex solutions to the equation /(a;) = that lie close to 
the line lies = D. First of all, such solutions must have e"'"^'^ close to e^"'^^, so we 
take u to be close to + 27rig/ wi, for an integer q. Then we write aq = p + x/ (27rz), 
for an integer p, which we will specify below, and uj = D + 2'Kiq/wi + A. With these 
substitutions, the equation fiuj) = transforms to l—e~^^^e~^^'^—e~'^'^^e~^e~'^'^^ = 

0. This equation defines A function of x. 

Lemma 6.2. Let Wi,W2 > and a = W2/W1 > 1; let D be such that ^-"^^0 j^^-w2D _ 

1, and let A = A(a;) be the function of x, defined implicitly by 

(6.6) ^-wW^-w^A ^ ^-nj2D^-x^-W2A ^ 

and A(0) = 0. Then A is analytic in x, in a disc of radius at least n around x = 0, 
with power series 

g-ui2D ^2g-uii_Dg-u!2-D 

A(x) = — :rT-^x+— 5 x^ + 0(x^), asx^O. 

^ ' f'{D) 2f'{Df ^ ' 



A PRIME ORBIT THEOREM FOR SELF-SIMILAR FLOWS 



21 



The coefficients of this power series are real. The coefficient of x is negative and that 
of is positive. 

Proof. Define y = y{x) by e~^^^y + Q-^^D^-Xya _ Then y{0) = 1 and wiA = 
— logy. Since y does not vanisli, it follows that if y{x) is analytic in a disc centered 
at X = 0, then A will be analytic in that same disc. Moreover, y is real- valued and 
positive when x is real, so A is real- valued as well when x is real. Further, y{x) is 
locally analytic in x, with derivative 

y [X) 



Hence there is a singularity at those values of x at which the denominator vanishes, 
which is at y = (a/(a — l))e"'^^ and = — a~"(Q; — Since this value is 

negative, the disc of convergence of the power series for y{x) is 

|x| < \—a logo; + (a — 1) log(a — 1) + 7ri| , 

which is a disc of radius at least vr. The first two terms of the power series for A{x) 
are now readily computed. □ 

Applying this, we find 
(6.7) uj = D + 2-ni^--—-x+^ % x^ + 0{x^), 

as X = 2'Ki{qa — p) ^0. We view this formula as expressing uj as an initial approxi- 
mation D + 27Tiq/wi, which is corrected by each term in the power series. The first 
corrective term is in the imaginary direction, as are all the odd ones, and the second 
corrective term, along with all the even ones, are in the real direction. The second 
term decreases the real part of u. 



Theorem 6.3. Let a be irrational and let p^ and he defined by (|6.3|) . Let q 

be a positive integer, and let q = Yl°^=k'^v<lv be the a-adic expansion of q, as in 
Lemma \6.]\ . Assume k > 2 or k = 1 and ai > 2, and put p = Yl^=k ^vPy. Then there 
exists a complex dimension of JF„ at 

q p-w2D 
uj = D + 2m 2Txi——{qa - p) 

f6 8) 

V ■ ^ ,2 ~wxD ~W2D 

The imaginary part of this complex dimension is approximately 2TTiq/wi, and its 
distance to the line Res = D is at least C/q'^^^^ where C = 2-K'^wle-^'"^^'"^^^ / f'{Df 
depends only on Wi and W2. 

Moreover, |Cw(s)| ^ q'k+2 around s = D+2niq/wi on the line Res = D, and \(n>{s)\ 
reaches a maximum of size C'{qa —p)~^, where C depends only on the weights Wi 
and W2. 
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Proof. By Lemma |6.1| , the quantity qa — p lies between {—^)^/(l'k+2 ^^'^ 

Under the given conditions on k, q'/^ > Qk ^ 2, hence x = 27ii{qa — p) is less than tt 

in absolute value. Then ( |6.7|) gives the value of u. 

Since the derivative of / is bounded on the line lies = D, this also implies that 
f{s) reaches a minimum of order {qa — pY on an interval around s = D + 27riq/wi 
on the line Res = D. It follows that |Cro(s)| ^ q'k+2 ^i^^ Res = D, with a 

maximum of order {qa — p)~'^. □ 

We obtain more precise information when q = qk- 

Theorem 6.4. For every k > {or k > 1 if ai = 1) , there exists a complex dimen- 
sion oj of J-'n, of the form 

(6.9, .^D, 2.^ - 2«(-l)' j;^^ - 2,=.; . O ) . 

as k oo. 

Moreover, \Cw{s)\ ^ q'k+i around s = D + 27iiqk/wi on the line Res = D, and 
\Cm{s)\ reaches a maximum of size C'q'^j^^, where C is as in Theorem 

Proof. Put p = pk- Then x = 27ri( — l)'^/g^_^^, which is less than vr in absolute value. 



The rest of the proof is the same as in the proof of Theorem |6.3| . □ 



Definition 6.5. A domain in the complex plane containing the line Res = D is a 

dimension free region for the flow J-'n, if the only pole of — Cro/Cw in that region is 
s = D. 

Corollary 6.6. Assume that the coefficients ao,ai, ... of a are bounded by M. Put 
B = Ti"^ e^^^^^'^'^'^^ / {2 f {D)^) . Then JF„, has a dimension free region of the form 

(6.10) {^ + ^^eC:a>D-^ 



oo 
fc=0 



The function —Cm/Cm satisfies hypotheses (Hi) and (H2) with k = 2. 

More generally, let a: M'*' [1, 00) be a function such that the coefficients {a^} 
of the continued fraction of a satisfy a^+i < a{qk) for every k > 0. Then J-'n, has a 
dimension free region of the form 

(6.11) {a + it eC: a > D 



t^a^{twi/{2TT)) 

If a grows at most polynomially , then —C'^/Qn satisfies hypotheses (Hi) and (H2) 
with K such that > t^a'^{twi/ {2-n)). 

Proof. This follows from Theorem |6.4| , if we note that for t = 2'n'qk/wi, we have 
q'k+i = c^k+iQk ^ '^a{qk)qk < ^a{qk)qk- So the complex dimension close to D + it is 
located at D -\- i{t -[- 0(g^"^\)) — (iyf/7r^)i?g^"Ij_^j^ + 0{qf7^^, where the orders denote 
real-valued functions. The real part of this complex dimension is less than D — 
Bt-'^a-\twi /{2n)). □ 

This has the following consequence for the Prime Orbit Theorem: 
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Theorem 6.7 (Prime Orbit Theorem with Error Term, for BernouUi flows). Let a 

W2/W1 have bounded coefficients in its continued fraction. Then 

(6.12) 'PM) = '4 + o(.'>('^^V"\. 



D \ \ \ogx 

as X 00. 

// a is 'polynomially approximable ', with coefficients in its continued fraction sat- 
isfying Ok+i < a{qk), for some increasing function a with a{x) = 0(a;'), as x ^ 00, 
then 

(6.13) Mx) = ^ + 0^x 

as X 00. 

The proof wiU be given in the next section, see Theorem |6.14 



Example 6.8. For the Golden flow, we have a = (p and wi = log 2 (see Exam- 
ple p.23| ). The continued fraction of is [1, 1, 1, . . .], hence = -Ffc+i, the {k + l)-st 
Fibonacci number, and = 0'^. Numerically, we find D ?a .7792119034 and the 
following approximation of A(x): 

- .47862 X + .08812 x^ + .00450 x^ - .00205 - .00039 x^ + .00004 + . . . . 

For every k, we find a complex dimension close to -D + 27rigfc/log2. For example, 
qg = 55, and we find a complex dimension at D — .00023 + 498. 58i. More generally, 
for numbers like g = 55 + 5 or g = 55 — 5 = 34 + 13 + 3, we find a complex dimension 
close to D + 27rig/log2, in this case respectively ai D — .023561 + 543. 63i and at 
D — .033919 + 453. 53z. In both these cases, the distance to the line Re s = D is 
comparable to the distance of the complex dimension for g = 5 to this line, which is 
located at Z) — . 028499 + 45. 05z. See Figure]^, where the markers are at the Fibonacci 
numbers. The pattern persists for other complex dimensions as well. Indeed, every 
complex dimension repeats itself according to the Fibonacci numbers. 

6.3. More than Two Generators. The following lemma replaces the continued 
fraction construction. 

Lemma 6.9. Let wi,W2, ■ ■ ■ ,wiy be weights such that at least one ratio Wj/wk is 
irrational. Then for every Q > 1, there exist integers q < Q^~^ and pj such that 
\qwj — PjW\\ < Wi/Q for j = 1, . . . , N. In particular, \qwj — PjWi\ < Wiq~^^^'^~^^ for 
j = l,...,N. 

Remark 6.10. Note that the condition implies that at least one ratio wj/wi is irra- 
tional. Also, \qwj — pjWi\ 7^ when wj/wi is irrational, so g cxd when Q ^ 00. 

The construction of such integers g and pj is much less explicit than for N = 2, since 
there does not exist a continued fraction algorithm for simultaneous approximation. Q 

^However, the L'^-algorithm can be used as a substitute for the continued fraction algorithm. 
We thank H. W. Lenstra, Jr. for guiding us to the following information: The L^-algorithm |LLL 



can be used to find fractions pj/q that approximate Wj/wi for j — 1, . . . , N . This algorithm works 
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The number Q plays the role of q'l^j^i in Theorem |6.4| above. In particular, if q is often 
much smaller than Q, then Wi, . . . , wjy is well approximable by rationals, and we find 
a small dimension free region. 

Again, we are looking for a solution of /(c<j) = close to s = D + 27riq/wi, where / 
is defined by ( |6.1|) . We write u = D + 27iiq/wi + A and wjq = wiPj+wiXj/ (27ri). For 
j = 1, we take pi = q and consequently xi = 0. In general, xj = 2m[qwj/wi — pj). 
Then f{uj) = is equivalent to 1 - e-'^^^e"^^"'"^^ = 0. 



The following lemma is the several variable analogue of Lemma |6.2| . However, 
in this case we do not know the radius of convergence with respect to each of the 
variables involved. 

Lemma 6.11. Let wi < W2 < ■ ■ ■ < wn, let D he such that XljLi ^""'■''^ ~ 1' "^^^ 
let A = A(x2, . . . , xn) be implicitly defined by 

N 

(6.14) J2 €-'"^^6-^^-'"^'^ = 1, 

i=i 

and Xi = 0. Then A is analytic in X2, ■ ■ ■ , x^, with power series 

(6.15) ' ' 

This power series has real coefficients. The terms of degree two form a positive 
definite quadratic form. 

Proof. The positive definiteness follows from the fact that the complex dimensions 



lie to the left of Res = D, see Theorem |2.19| . It can also be verified directly. □ 
Applying this, we find 

UJ = D + 2TXi > ; zXj-\ > ; X? 

j^J'{D) ' 2^ f'iD) ^ 
(6.16) ^ ' ' 

where Xj = 2TTi{qWj/wi —Pj). Again, this formula expresses u as an initial approx- 
imation D + 2'n'iq/wi, which is corrected by each term in the power series. The 
corrective terms of degree one are again in the imaginary direction, as are all the odd 
degree ones, and the corrective terms of degree two, along with all the even ones, are 
in the real direction. The degree two terms decrease the real part of u. 



in polynomial time, like the continued fraction algorithm, but it does not give the best possible 
value for q (given a certain error of approximation). The problem of finding the best value for q is 



NP-complete |Lag|. (See also [GL£].) 
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Theorem 6.12. Let N > 2 and let Wi, . . . ,W]sr be weights. Let Q and q be as in 

Lemma \6.9{ . Then JF„ has a complex dimension close to D + 27iiq/wi at a distance 
of at most 0{Q~^) from the line lies = D, as Q ^ oo. The function \Cm/Cro\ reaches 
a maximum of order Q^. 

Proof. Again, the numbers Xj are purely imaginary, so the corrective terms of degree 1 
(and of every odd degree) give a correction in the imaginary direction, and only the 
corrective terms of even degree will give a correction in the real direction. Since 
\xj\ < 2tt/Q, the theorem follows. □ 

Corollary 6.13. The best dimension free region that J^^ can have is of size 

(6.17) {a + 2t: a >D-0 (r 2/(^-1))}. 

The implied constant depends only on wi, . . . ,wn- 
If Wi, . . . , wn is 'a-approximable' , then the dimension free region has the form 

(6.18) {a + it: a>D~0 {a-^{wit / {2Ti))t-^/^''-^^) ] , 

where a: [l,oo) is an increasing function such that for every integer q > 1, 

\qWj — PjWil > {wi/ a(q))q~^^^^~^^ for j = 1, . . . , N. 

This has the following consequence for the Prime Orbit Theorem: 

Theorem 6.14 (Prime Orbit Theorem with Error Term). Suppose Wi, . . . ,wn are 

badly approximable, in the sense that \qwj — PjWi\ ^ for j = 1, . . . , N 

and every q > 1. Then 



^a^n\ i ( \ ^ n ( D /log log X 

(6.19) ^m{x) = — + 0\x 



N-l ■ 
4 



D \ \ logx 
as X ^ oo. 

// Wi,...,WAr is 'polynomially approximable', in the sense that \qwj — PjWi\ > 
{wi / a{q))q~^^^^~^^ for j = 1, . . . , N and every q > 1, for some increasing function a 
on [1, oo) such that a{x) = O(x') as x —>■ oo, then 

, f ^ D /loglogx\ *'(iv-i)+4 

(6.20) Mx) = ^ + 0\x y-^^^] 

as X — > oo. 



Proof of Theorems \6. Tj and \6.14 - We apply the pointwise explicit formula at level k 



2 (see Theorem p.l|) to obtain 

The error term is estimated by i?[^l(x) = 0{x^^^^'') for some positive c. We will 
estimate the sum by an argument which is classical in the theory of the Riemann 
zeta function and the Prime Number Theorem, under the assumptions that the u 
have a linear density, and that every u = a + it satisfies a < D — Ct^^ for some 
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positive number p. We then obtain Theorem |6.14| by taking p = 2/{N — 1) + 21, and 
Theorem |6.7| corresponds to the case when N = 2. 

The sum Jl^+i) absolutely convergent. We split this sum into the parts with 

I Imto"! > T and with | lmu\ < T. Put A = \u!{u! + From the fact that 

the complex dimensions have a linear density, it follows that there exists a constant 



B is such that E|imo.|>T 
^^n+i-CT-^ + 5x^+VT. For T 



ti;(ti; + l)| < B /T for every T. Then ^ 
(pClogx/ log logo;) we find 



< 



u{iu + 1) 



O 



D+l /^loglogX 

logx 



We then apply a Tauberian argument to deduce a similar error estimate for 



see M, p. 64]. Let h = x (log logo;/ logo;) ^/'^^^^ Thus 



h 



x+h 



1pm{t) dt 



h 



[2], 



Now 



(x+h) 



D+l_ 



hD{D+i) ~ = x^/D + 0{x^-^h) = x^ / D + x^O{iloglogx/ logx 
ther, 0{x^^^{loglogx/ logxY^f' /h) = a;^0((logloga;/logx)^/(^^)). 



|l/(2p)^ 



Fur- 
□ 



Remark 6.15. Note that by using the Tauberian argument, we lose a factor two in 
the exponent. Indeed, the estimate 

JV-l 



y ^- + R{x) = x^C-^f^) 



2i(]V-l)+2 



holds distributionally. 

Remark 6.16. If a{q) grows more than polynomially, we obtain a bound of the form 
x'^/a™^(logx) for the error in the Prime Orbit Theorem, where a™^ is the inverse 
function of a. 

6.4. Conclusion. If / > in the exponent (iV — l)/(4/(A^— 1)+4) of log x in the error 
term of Theorems and |6.14|, then the error term is independent of A^, essentially 



of order x'^ (log x) "^^^ (ignoring the factor of log logx). Thus, if the weights are well 
approximable, the error term is never better than x^ divided by a fixed power of the 
logarithm of x. On the other hand, when / = 0, that is, roughly speaking, when the 
weights are never close to rational numbers, then the error term is essentially of order 
x'^(logx)~*^^~^'*/''. Hence, the larger iV, the smaller the error term in that case. 

We may compare this, somewhat superficially in view of the Riemann Hypothesis, 
with the situation of the Riemann zeta function. In view of Example |3.6| , the weights 
are Wp = tt)j(p) = logp, for each prime number p, and there are infinitely many of 
them. Since it is expected that {logp}p:prime is badly approximable, one expects an er- 
ror term of order "x'^(logx)~°°" . Indeed, in (|3.22| ), we have e~'^^^°^^ = O ((logx)~^) 
for every > 0. The corresponding pole free region has width A/ logt at height t 
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(see [I, Theorem 19]), which is "t This lends credibihty to the conjecture that 

{logp}p:primc IS badly approximable by rational numbers. 

Acknowledgement. The authors wish to thank Gabor Elek for helpful conver- 
sations about dynamical zeta functions. 
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